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The interplay between crystallinity and super- 
fluidity is of great fundamental and technolog- 
ical interest in condensed matter settings. In 
particular, electronic quantum liquid crystallinity 
arises in the non-Fermi liquid, pseudogap regime 
neighboring a cuprate's unconventional supercon- 
ducting phase^. While the techniques of ultra- 
cold atomic physics and quantum optics have 
enabled explorations of the strongly correlated, 
many-body physics inherent in, e.g., the Hubbard 
modeP, lacking has been the ability to create a 
quantum degenerate Fermi gas with interparti- 
cle interactions — such as the strong dipole-dipole 
interactiorF' capable of inducing analogs to elec- 
tronic quantum liquid crystals. We report the 
flrst quantum degenerate dipolar Fermi gas, the 
realization of which opens a new frontier for ex- 
ploring strongly correlated physics and, in par- 
ticular, the quantum melting of smectics in the 
pristine environment provided by the ultracold 
atomic physics setting"*. A quantum degenerate 
Fermi gas of the most magnetic atom ^^^Dy is 
produced by laser cooling to 10 /iK before sym- 
pathetically cooling with ultracold, bosonic *^^Dy. 
The temperature of the spin-polarized ^^^Dy is a 
factor T/Tp — 0.2 below the Fermi temperature 
Tp = 300 nK. The co-trapped ^"^^Dy concomitantly 
cools to approximately Tc for Bose-Einstein con- 
densation, thus realizing a novel, nearly quantum 
degenerate dipolar Bose-Fermi gas mixture. 

Quantum soft phases are states of quantum matter 
intermediate between canonical states of order and dis- 
order, and may be considered the counterparts of liq- 
uid crystalline and glassy states in classical (soft) con- 
densed matter physics. Such phases tend to arise under 
competition between short and long-range interactions 
and often result in the non-Fermi liquid, strongly corre- 
lated behavior manifest in some of the most interesting 
electronic materials of late: high- Tc cuprate supercon- 
ductors, strontium ruthenates, 2D electron gases, and 
iron-based superconductor^. Recent theory suggests the 
long-range, anisotropic dipole-dipole interaction (DDI) 
among atoms in a degenerate Fermi gas, confined in an 
harmonic trap or optical lattice, may also induce transi- 
tions to states beyond the now-familiar insulating, metal- 
lic, and superfluid. Namely, phases that break rotational, 
translational, or point group symmetries may emerge in 
a manner akin to those found in classical liquid crystals. 



e.g., the nematic and smecticP. 

Degenerate gases of highly magnetic fermionic atoms, 
such as ^^^Dy, may shed light on QLC physics without 
unwanted solid state material complexity, disorder, and 
dynamical lattice distortions. Uniaxial (meta-nematicj^ 
and biaxial nematicf^ distortions of the Fermi surface of 
a harmonically trapped gas in the presence of a polariz- 
ing field may be observable as well as meta-nematic and 
smectic phases in 2D anisotropic optical latticed^H^. An 
exciting prospect lies in the possibility of achieving spon- 
taneous mag netiza tion in dipolar systems coupled with 
nematic orde^^^^^. Additionally, DDI-induced pairing 
of fermions may lead to supersoliditjB^ and bond order 
solidiP. 

However, obtaining a quantum degenerate dipolar 
Fermi gas has been a difficult, unrealized experimental 
challenge. The highly magnetic fermionic atoms ^'^Cr (6 
Bohr magnetons ^j, b) an d *^^Er (7 /is) have yet to be 
cooled below 10 /xKEIini. The fermionic polar molecule 
"'"K^^Rb (0.57 Debye) has been cooled to near degener- 
acy {T /Tp — 1.4|i3 and loaded into a long-lived lattice 
while partially polarized (0.2 D^^^l^ b^^t complexities aris- 
ing from ultracold chemistry have hampered additional 
evaporative coolinj^^. In contrast, magnetic fermionic 
atoms do not undergo chemical reactions and are im- 
mune to inelastic dip olar c ollisions when spin polarized 
in high magnetic fieidplani. 

The strong, character of the DDI arises in ground 
state polar molecules though a polarizing electric field 
that mixes opposite parity states. This electric field 
breaks rotational symmetry; consequently, observing the 
full range of true (non-meta) quantum nematic and fluc- 
tuating smectic phases, and their often unusual topolog- 
ical defects, is not possible in systems of fermionic polar 
molecules, especially in three dimensions. By contrast, 
highly magnetic atoms exhibit the DDI interaction even 
in the absence of a polarizing field. Moreover, the mag- 
netic DDI can be tuned from positive to negativ^^Il^ which 
may be important for simulating dense nuclear matter. 

Dysprosium's isotopes ^^^Dy and *^^Dy are the most 
magnetic fermionic atoms. With a dipole moment of /i = 
10 /iB, *^^Dy provides a DDI length Zddi = iLQ^'^m/ Aiih^ 
that is factors of [400, 8, 2] larger than that of ["OR, ^^Cr, 
i67Er]. With respect to fully saturated ^OR^^Rb (0.57 D), 
^^-'^Dy is 30x less dipolar for equal densities, but within 
a factor of 2 if confined in a lattice of less than half the 
periodicity. Lattices of wavelength 400-500 nm may be 
possible with Dy, whereas for molecules, photon scat- 
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FIG. 1. Dy transitions and laser cooling and trapping scheme, a, Electronic energy level structure for bosonic Dy 
(nuclear spin 7 = 0) and fermionic ^^^Dy (/ — 5/2), including laser cooling and trapping transitions, b, Additional ground 
and excited-state hyperfine structure exists for ^'^^Dy (F — I + J, where J = 8 is the total electronic angular momentum and 
primes denote the excited states). Shown is the 32-MHz-wide transition at 412 nm used for the transverse cooling, Zeeman 
slower, capture MOT, and imaging beams, c. Blue-detuned, narrow-line (1.8 kHz-wide) MOT cooling transition at 741 nm. 
d. Spectra of 421-nm '^^^Dy and ^^^Dy (hyperfine) transitions including relative detunings of each MOT and Zeeman slower 
(ZS) laser, e. Transition and detuning spectra for MOT on 741-nm line, f , Sketch of dual species crossed optical dipole trap 
(ODT) aspect ratio along with magnetic field B and gravity g orientation. The imaging (421-nm) beam and the orthogonal 
ODT (1064-nm) beams are in the p- plane. 



tering from rovibronic states at these wavelengths may 
reduce gas lifetimes. Indeed, ^^^Dy confined in a 450- 
nm lattice would provide a DDI strength more than 3x 
larger than ""^K^^Rb confined in a 1064-nm lattice, if the 
electric dipole moment is unsaturated (0.2 D) to maintain 
coUisional stability We estimate that the DDI strength 
of Dy confined in short-wavelength lattices would be suf- 
ficient to observe some of t he exotic many-body physics 
recently proposed^'' 2223]^ 

Until recently, the laser cooling of Dy posed an insur- 
mountable challenge due to its complex internal struc- 
ture and the limited practicability of building repumping 
lasers: ^140 metastable states exist between the ground 
state and the broadest laser cooling transition at 421 nm 
(see Fig. [Ti). Moreover, an open /-shell submerged un- 
derneath closed s-shells, combined with a large magnetic 
moment and electrostatic anisotropy from the L = 6 
orbital angular momentum, pose challenges to molecu- 
lar coupled-channel calculations^ which could otherwise 
guide early experiments. Even more daunting is the ad- 
ditional hyperfine structure of the fermions, shown in 
Fig. [TJj-c, which splits each level into six. 

Despite this vast energy level state-space, a repumper- 



less magneto-optical (MOT) technique was able to indi- 
vidually laser cool and trap of all five high-abundance 
isotopes (three bosons, two fermions) to 1 mK us- 
ing a single laseil^S (see Methods for details). More- 
over, the complex homonuclear molecular potentials — 
involving 153 Born-Oppenheimer surfaces — and the as- 
sociated multitude of scattering lengths did not inhibit 
the efficient Bose-condensation of spin-polarized ^^''Dy 
through forced evaporative cooling^^. 

However, quantum degeneracy of identical Fermi gases 
is often more difficult to achieve than Bose-condensation 
because s-wave collisions are forbidden due to the re- 
quirement that the total wave function for two identi- 
cal fermions be anti-symmetric with respect to particle 
exchange. Rethermalization from elastic collisions cease 
below the threshold for p-wave collisions (at typically 10- 
100 ^J-K)^ and efficient evaporative cooling can no longer 
be maintained. Co-trapping mixtures of particles — either 
as different spin states of the same atom or as mixtures 
of isotopes or elements — reintroduces s-wave collisions, 
providing a finite elastic cross-section for scattering even 
down to low temperatures. The mixture needs to be sta- 
ble against inelastic collisions which could add heat or 



3 





A 


< 


b 

O , 



FIG. 2. Images of the Dy degenerate Fermi gas. a, 

Single sliot time-of-fliglit absorption image at t — 6 ms. b, 
Average of six images. Density integrations versus p (c) and 
z (d). Tlie green curve is a gaussian fit to tiie data's wings 
(radius a — 20 ^m), while the red curve is a fit to a Thomas- 
Fermi distribution. Data are consistent with a Thomas- Fermi 
distribution of T/Tf = 0.21(5). The Fermi velocity and tem- 
perature are 5.6(2) mm/s and 306(20) nK, respectively, and 
the gas temperature is 64(16) nK. The degenerate Fermi gas 
contains 6.0(6) x 10^ atoms at peak density 4(1) x lO'^^ cm"^. 



induce trap loss. Evaporating ^^-"^Dy in a mixed state 
of two spins, as proved efficient for 40j^|26l would lead to 
large dipolar relaxation-induced heating even in the pres- 
ence of small, mG-level fields because the inelastic, sin- 
gle spin- flip cross-section cti — aC,{kf/ki) scales strongly 
with dipole momentP^ 

^"15 ^^^2 4nh^ ) 

where Fi is the spin of the atom whose spin flips (-Fi = F2 
for identical particles), gi are (7- factors for atom i, m is 
mass, and ki and kf are the initial and final momenta. 
For leiDy (i62Dy), F ^21/2 {F = J ^ %) and gp {gj) 
= 0.95 (1.24). The function C,{kf/ki) = [l + eh{kf/h% 
where e = ±1, and h{x) is defined in Methods, accounts 
for the contributions of even, odd, or all partial waves to 
the scattering process. 

We choose, therefore, to seek a degenerate dipolar 
Fermi gas with Dy by sympathetically cooling ^^^Dy with 
the boson ^^^Dy while both are spin-polarized in their 
strong-magnetic-field seeking ground states: \F,mp) — 
|21/2,-21/2) for ^^iDy and \J,mj) ^ |8, -8) for ^^^^Dy. 
See Fig. [T^-c for energy level schemes. Preparation of 
this ultracold Bose-Fermi mixture — the first such mix- 
ture for strongly dipolar species — builds on our single- 
species techniqu^i^ for Bose-condensing ^^''Dy and relies 



on the laser cooling and trapping of two isotopes before 
loading both into an optical dipole trap (ODT) for forced 
evaporative cooling. We sketch here the experimental 
procedure; further details are provided in the Methods. 

Isotopes ^^^Dy and ^^^Dy are collected sequentially in 
a repumperless MOT operating on the 421-nm transi- 
tion25 (Fig. [l}i), with final MOT populations of = 
2 X 10'' and 4 x lO'', respective ly. N ext, simultaneous 
narrow-line, blue-detuned MOT^i^ias) (,qq^ both isotopes 
to 10 /xK via the 741-nm transition (Fig. [l^) for 5 s to 
allow any remaining metastable atoms to decay to the 
ground state . The blue-detuned MOTs also serve to spin 
polarizd^^l^ both isotopes to their maximally high-field- 
seeking (metastable) states mp = -\-F (nij — +J) for 
leiDy (i62Dy). 

The blue-detuned MOTs of the two isotopes can be 
spatially separated due to the d ependence of the MOTs' 
positions on laser detunin ^^^ l ^^ l This allows the isotopes 
to be sequentially loaded into the 1064-nm ODTl in 
Fig.fTf, which is aligned above the ^^^Dy MOT but below 
the ^^By MOT. First ^^"^By and then i^^Dy is loaded 
into ODTl by shifting the quadrupole center with a verti- 
cal bias field. All 741-nm light is extinguished before the 
spin of both isotopes are rotated via radiofrequency (RF) 
adiabatic rapid passage (ARP) into their absolute ground 
states mp = -F (rrij = -J) for ^^^Dy {^^^Dy). The 
ODTl populations of ^^^^Dy and ^^^Dy are both initially 
1 X 10^ before plain evaporation cools the gases to 1-2 /iK 
within 1 s. A 0.9 G field is applied close to the trap axis 
of symmetry z throughout plain and forced evaporation. 
This provides a Am = 1 Zeeman shift equivalent to 50 
(70) /iK for ^^^Dy (^^^Dy). Because this is much larger 
than the temperatures of the gases, the field serves to 
maintain spin polarization while stabilizing the strongly 
dipolar ^^^Dy Bose gas against collapse as its phase-space 
density increases'^'^. 

Magnetic Stern-Gerlach measurements and observa- 
tions of fluorescence versus polarization are consistent 
with an RF ARP sequence that achieves a high degree 
of spin purity for each isotope. Remnant population in 
metastable Zeeman substates quickly decays to the abso- 
lute ground state via dipolar relaxation regardless of col- 
lision partner at a rate of F cx ainv = 1-10 s~^, where n 
is the atomic density and v is the relative velocity during 
the plain evaporation stage. (Since gpF — gjJ, colli- 
sions between Bose-Bose, Bose-Fermi, and Fermi-Fermi 
pairs result in F's of similar magnitude as long as ki ~ k ^ 
since h{x^l)^Q. This condition is fulfilled during plain 
evaporation due to a low ratio of Zeeman-to-kinetic en- 
ergy. For example, inelastic dipolar ^^^Dy-^^^Dy colli- 
sions (e = —1) proceed at rate F = 1-5 s~^ even in the 
absence of ^^^Dy.) Thus, a (two-body) coUisionally sta- 
ble mixture of identical bosons and identical fcrmions is 
prepared within the 1 s between spin rotation and forced 
evaporation. 

Subsequently crossing ODTl with 0DT2 forms 
an oblate trap with frequencies [/xj/yj/z] = 
[500, 580, 1800] Hz. Ramping down the optical power 
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FIG. 3. Temperature and population of ^'^^Dy and ^"^Dy mixture, a, Three measures (see text) of ^'^^Dy p temperature- 
to-Fermi temperature Tf as trap population decreases due to forced evaporation of the spin-polarized Bose-Fermi mixture, 
b, The dipolar Bose-Fermi mixture remains in thermal contact throughout the evaporation sequence as measured by fits of 
TOF densities to gaussian or TF distributions. The orange dashed line demarcates the boundary below which the temperature 
of ^^^Dy is lower than Tf- Likewise, the purple dashed line demarcates the temperature below which Tc for ^^^Dy Bose 



degeneracy would be reached given the trap frequencies and population at 
Bose-Fermi mixture versus evaporation time. 



19 s. c, Trap populations of the spin-polarized 



lowers the trap depth and evaporates spin-polarized 
^^^Dy to quantum degeneracy; with a 19 s evaporation, 
the final trap has frequencies [180, 200, 720] Hz and 
id/2'!T is defined as their geometric mean. Figure [2] shows 
the density profile of ultracold ^^^Dy, which is more con- 
sistent with a Thomas-Fermi (TF) distribution arising 
from Fermi-Dirac statistics than a gaussian arising from 
a classical, Maxwell-Boltzmann (MB) distribution. 

Figure [3] shows a collection of such fits as a function of 
trap population and evaporation time. The T/Tp data 
are extracted from density profiles using three methods: 
1) The high momentum wings are fit to a MB distribution 
with ksTp = hui{6Ny^^ extracted from measured trap 
parameters and population (Fig. |3j:); 2) similar to (1), 
but using a TF distribution; 3) fitting the fugacity to 
directly extract T/Tp. The last method is known to be 
inaccurate at higher temperatures, while the gaussian fit 
tends to overestimate the temperature below T/Tp = 
LfSslThe evaporation does not yet seem to reach a plateau 
in cooling efficiency at 19 s; poor imaging signal-to-noise 
hampers measurements at longer evaporation times. 

Data in Fig. show that thermal equilibrium be- 
tween the bosons and fermions is maintained through- 
out the evaporation, and Bose-condensation of ^^^Dy 
within the mixture is nearly reached for an evaporation 
of 19 s. We estimate the corresponding critical temper- 
ature Tc « 40 nK of co-trapped -'^^^Dy by scaling the 
measured ~ 120 nK of singly trapped ^^^Dy with the 
cube root of their relative trap populations. (^^^Dy has 



been Bose-condensed in the absence of ^^^Dy; manuscript 
in preparation.) The nearly doubly degenerate dipolar 
Bose-Fermi mixture may lead to interesting dipolar and 
many-body physics once cooling efficiency improves. 

While interacting BECs invert their anisotropic aspect 
ratio upon time-of-flight expansion, anisotropic degen- 
erate Fermi gases tend to a spherical shape. As the 
DDI strength increases, the degenerate Fermi gas will 
expand into a prolate ellipsoid oriented along the mag- 
netization direction regardless of trap aspect ratio. Fur- 
thermore, the gas may become unstable when the quan- 
tity edd = ^iQ^i'^{m?C) /l&n'^hl'Y''^ > l.I^ At the lowest 
attained T/Tp, edd — 0.2 for ^^^Dy, and the ratio is 
Iddi/If is 0.05. This DDI strength should lead to Fermi 
surface distortions (as yet unmeasured) at the percent 
levepfi-. Both ratios could be enhanced ^3x by increas- 
ing trap frequency using a more spherical confinement — 
while maintaining stability of the dipolar Bose gas — and 
by increasing trap population. Additionally, the effi- 
ciency of evaporation — from gaussian fits to the data, 
7 = 3 '^^^"j^^^'' ~ 2.3 — may be improved by optimally 
tuning the magnetic bias field near one of the presum- 
ably large number of Feshbach resonances. 

Surprisingly, we achieve the forced evaporative cooling 
of spin-polarized ^^^Dy without ^^^By to T/Tp = 0.7 at 
Tp — 500 nK. As mentioned above, achieving quantum 
degeneracy with spin-polarized identical fermions alone is 
usually not possible due to suppression of elastic scatter- 
ing below the p-wave threshold 50 fiK.^'^ That such a low 
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temperature ratio is achieved may be a novel consequence 
of the highly dipolar nature of this gas: namely, that a 
significant elastic cross-section persists to low tempera- 
tures due to the as yet unobserved phenomenon of univer- 
sal dipolar scatterin^SI. The associated scattering rate 
is expected to scale as m'^/^/i^ regardless of the details 
of the short-range molecular potential. The predicted 
fermionic Dy universal dipolar cross-section, 7.2 x 10~^^ 
cm^, is nearly equal to ®^Rb's s-wave cross-section and 
could provide sufficient rethermalization for the evapo- 
rative cooling we observe. While future measurements 
will quantify spin purity, Fermi statistics does not in- 
hibit identical fermions from spin purifying via dipolar 
relaxation once the RF ARP sequence populates the ma- 
jority of the atoms in the absolute ground state. We will 
present these data along with supporting measurements 
of coUisional cross-sections and scattering lengths in a 
future manuscript. 

The efficacy of experimental proposals^ for studying 
the quantum melting of QLCs — important for better un- 
derstanding the relationship of incipient electronic stripe 
order to unconventional superconductivity — may now be 
investigated using degenerate ^^^Dy in optical lattices 
wherein Iddi/If can be greatly enhancetP^. Looking 
beyond QLC physics, the large spin F — 21/2 of the 
novel degenerate dipolar Fermi gas presented here opens 
avenues to explore exotic spinor physics as well as physics 
associated with strong spin-orbit coupling. 

I. METHODS 

Repumperless MOT. The 421-nm Zeeman slower laser 
is detuned -650 MHz from resonance, and a 421-nm MOT 
collects 2 X 10^ ^^^Dy atoms in 4 s with the aid of a 
pre-slower transverse cooling stage. The MOT operates 
without repumpers despite many decay channels because 
the highly magnetic, metastable Dy remains trapped in the 
magnetic quadruple trap for a sufficiently long time to decay 
to the ground statd^. All beams derived from the 421-nm 
laser are then shifted down in frequency by 1.18 GHz and a 
^'^^Dy MOT collects 4 x lO'' atoms in 100 ms. This scheme 
balances the need for preservation of ^''^Dy population 
against coUisional decajf=^ with the trapping of sufficient 
^^^Dy for subsequent sympathetic cooling. Reversing this 
procedure results in fewer atoms due to heating of ^^'^Dy in 
the magnetic trap by light scattered from the ^^^Dy Zeeman 
slower beam. 

Narrow-line MOT. The narrow-line MOT lasers are blue- 
detuned 0.6 MHz from the 2-kHz wide 741-nm transition. 
The MOTs form below the quadrupole center such that 
the atomic transition is Zeeman-shifted to the blue of the 



laser. The MOT positions are determine d by t he balance of 
optical, magnetic, and grativational forcej^^"^, and because 
the laser detunings for ^^^Dy and ^^^Dy can be adjusted 
independently, the two clouds can be easily separated. 

ODT, RF ARP and crossed ODT. ODTl has a waist 
of 30 (60) urn and ODT2 has a waist of 20 (60) fim in z 
(p). ODTl is ramped on in 100 ms after 5 s of narrow-line 
coohng. After loading ODTl— 120 (800) ms for ^"^^Dy 
(^"^Dy) — a 4.3 G field is applied while a 20-ms RF sequence 
flips the spin of both isotopes to the absolute ground state. 
Dipolar relaxation further spin purifies the gas, and within 

the cross section ai, h{x) = - I ''^(ifj] log ((TT^) ^^^^ 
ReiW\). After spin purification, 0DT2 is turned on, forming 
a crossed ODT with depth 300 fj,K. ODTl and 0DT2 have 
initial powers of 18 W and 12 W. The background limited 
lifetime exceeds 180 s. 

Forced evaporation scheme. The two beams of the 
crossed ODT are ramped down according to the functional 
form P{t) — -Po/(l + */''")'' using experimentally determined 
parameters r = 6 s and /3 = 1.5 for both beams. Forced 
evaporation with 0DT2 begins 5 s after ODTl. The 
beam powers are held at the final values for 400 ms after 
evaporation to allow for equilibration. 

Temperature fitting. The TOP images are fit to both gaus- 
sian and TF profiles, the latter of the form^'^ 

.Z^J^H"^^ "^]^" (1) 

with fitting parameters A, yo, xo, ay, a^, and c. The 
fitting is scaled by the constant Li2(— C) such that param- 
eter A corresponds exactly to the image peak OD. Because 
the high-velocity component of the cloud is less sensitive to 
signatures of quantum degeneracy, only the wings of the ex- 
panded cloud are used for the gaussian TOF fits. From the 
TF fits, the temperature is determined from size of the cloud, 
o"TF.i = {kbT/muJi)[l + {LOit)'^], where uji is the trap frequency 
in 2 or p and t is the expansion time. The cloud fugacity ( 
is also determined by the TF fits and provides a direct mea- 
sure oiT/Tp = [-6Li3(-C)]"^/^, where Lis is the third-order 
polylogarithm function. 
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